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We measure the cross-correlation of Atacama cosmology telescope cosmic microwave background 
(CMB) lensing convergence maps with quasar maps made from the Sloan Digital Sky Survey DR8 SDSS- 
XDQSO photometric catalog. The CMB lensing quasar cross-power spectrum is detected for the first time 
at a significance of 3.8er, which directly confirms that the quasar distribution traces the mass distribution 
at high redshifts z > 1. Our detection passes a number of null tests and systematic checks. Using this 
cross-power spectrum, we measure the amplitude of the linear quasar bias assuming a template for its 
redshift dependence, and find the amplitude to be consistent with an earlier measurement from clustering; 
at redshift z ~ 1.4, the peak of the distribution of quasars in our maps, our measurement corresponds to a 
bias of b = 2.5 ± 0.6. With the signal-to-noise ratio on CMB lensing measurements likely to improve by 
an order of magnitude over the next few years, our results demonstrate the potential of CMB lensing cross- 
correlations to probe astrophysics at high redshifts. 
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I. INTRODUCTION 

As the photons of the cosmic microwave background 
(CMB) travel through the Universe, they are gravitation- 
ally deflected by the web of matter through which they 
pass. In the CMB sky we observe today, these deflections 
are imprinted as arc-minute-scale distortions of small- 
scale temperature fluctuations [1,2]. The microwave back- 
ground thus contains information not only about the 
primordial Universe at redshift ~ 1 100, but also about 
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the matter density fluctuations in the lower-redshift 
Universe. 

The lensing deflection field d(n) points from the direc- 
tion n in which a CMB photon was received to the direc- 
tion from which it was emitted. This deflection field can be 
determined from the measured CMB because lensing 
changes the statistics of small-scale unlensed CMB fluctu- 
ations in a characteristic way, introducing correlations 
between different Fourier modes. By measuring correla- 
tions between pairs of Fourier modes that would be 
uncorrelated in the absence of lensing, one can estimate 
d [3] and hence the lensing convergence k = — V • d/2 
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(a useful quantity because it is a direct measure of the 
projected matter density, see Eq. (1)). Using a quadratic 
estimator, lensing was first measured in cross-correlation 
with radio sources and galaxies by Refs. [4,5] (using 
WMAP data) and in autocorrelation by Ref. [6] (using 
Atacama cosmology telescope data). More accurate mea- 
surements of both the lensing power spectrum and the 
lensing galaxy cross-correlation were recently reported 
by the South Pole Telescope [7,8]. 

CMB lensing measurements are a powerful cosmologi- 
cal probe [9] because they are sensitive to both the growth 
of density fluctuations and the geometry and size of the 
Universe, yet are relatively insensitive to both instrumental 
and astrophysical systematic errors [6,7]. Lensing mea- 
surements can constrain the properties of dark energy 
[10], the amplitude of density fluctuations and the masses 
of neutrinos [11]. They can also constrain the properties of 
biased tracers of the matter distribution. The focus of this 
paper is the cross-correlation of CMB lensing maps with 
one such tracer — quasars. 

Quasars are among the most luminous objects in the 
Universe. Their immense luminosity is believed to be pow- 
ered by accreting supermassive black holes [12,13] which 
reside at the center of almost every massive galaxy [14]. As 
the activity of quasars and the rate of star formation appear 
to be linked [15], they are a crucial element in our present 
understanding of galaxy evolution. Measurements of the 
relation between dark matter and the distribution of quasars 
can inform us about quasar properties such as the masses of 
the dark matter halos that host the quasars, the scatter in the 
quasar-halo mass relation and the quasar duty cycle (see 
e.g., Ref. [16]). Such measurements of quasar properties 
will, in turn, improve our understanding of structure for- 
mation and galaxy evolution. 

Both the number density of quasars and the strength of 
CMB lensing in a certain direction depend on the projected 
dark matter density in this direction, and quasars are most 
common at the redshifts that produce the largest lensing 
deflections. This implies that the CMB lensing and quasar 
fields should be strongly correlated [17]. Measuring the 
cross-power spectrum and comparing it to theoretical cal- 
culations, we can determine the proportionality factor 
which relates a fluctuation in matter density to a fluctuation 
in quasar number density. This proportionality factor is 
known as the quasar linear bias b (defined as b = 8 / 8, 
where 8 q and 8 are the fractional spatial overdensities of 
quasars and matter, respectively). 

In this work we present the first measurement of the 
CMB lensing quasar cross-power spectrum, and use it to 
derive a constraint on the quasar bias. The paper is organ- 
ized as follows: Sec. II presents the theoretical background 
underlying the CMB lensing quasar cross-correlation. 
Section III explains how the lensing and quasar maps 
used in our analysis are constructed, and describes the 
resulting cross-spectrum measurement. The constraint on 


quasar linear bias we obtain from the cross-power spectrum 
is presented in Sec. IV. Null tests and systematic checks are 
discussed in Sec. V. All calculations assume a ACDM 
cosmology with WMAP5 parameters [18] and cr 8 = 0.819. 

II. THEORETICAL BACKGROUND 

Cosmological weak lensing effects can be described 
using the convergence field k, which is equal to a weighted 
projection of the matter overdensity 8 [1]: 

C Z-LS 

k( n ) = / ' dzW K (z)8(r](z)n, z). ( 1 ) 

Jo 

The relevant kernel (assuming a flat universe) is 
3 (ti^s — 77 ( 7 )) 
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where 17 ( 2 ) is the comoving distance to redshift z, n is a 
direction on the sky, r/ LS is the comoving distance to the 
last scattering surface, zls i s the redshift of the last scat- 
tering surface, H(z) is the Hubble parameter, and fl 0 and 
// 0 represent the present values of the matter density 
parameter and the Hubble parameter, respectively. 

The fractional overdensity of quasars in a direction n is 
given by </(n), which — assuming a linear bias relation be- 
tween the distribution of matter and quasars — is given by 

q( n) = / ' dzW q (z)S(r)(z)n, z ), (3) 

Jo 

where the kernel is 
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and 

, dN 

g(z) = viz) J , S dz!{ 1 - viz)/ viz 1 ) ) |-y^ rfjyj ( 5 ) 

(see Ref. [17]). Here b is the linear bias, dN / dz is the redshift 
distribution of quasars, and .v is the variation of the 
number counts of quasars N(<m ) with the limiting 
magnitude m at the faint limit of the quasar catalog, 
.v = d\og A(] N / dm. The second term in Eq. (4) is the magni- 
fication bias, the change in source density in a certain direc- 
tion induced by lensing magnification. For the quasar catalog 
used in this work (which has s ~ 0.043), this term is sig- 
nificantly smaller than the first term (~ 15% of its magni- 
tude) and is negative. 

Combining Eqs. (1) and (3) gives the expected 
lensing quasar cross-power spectrum in the Limber 
approximation [19]: 

c ? = f W K (z)W‘'(z)P(k = i/viz), z), ( 6 ) 
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where P(k, z) is the matter power spectrum at wavenum- 
ber k and redshift z. 

III. CROSS-CORRELATING CMB LENSING 
AND QUASARS 

A. The ACT CMB lensing convergence maps 

The lensing convergence fields used in our analysis are 
reconstructed from CMB temperature maps made by the 
Atacama cosmology telescope (ACT) [20-22], a 6-meter 
telescope operating in the Atacama desert in Chile. These 
CMB temperature maps are obtained from observations 
made during 2008-2010 in the 148 GHz frequency band 
and calibrated as in Ref. [23]. The maps consist of six 
patches, each of size 3X18 degrees, in the Sloan Digital 
Sky Survey (SDSS) Stripe 82 region [24], with a map- 
averaged white noise level of 21 /x K arc min. Radio and 
IR point sources as well as Sunyaev-Zel’dovich (SZ) clus- 
ters detected with a matched filter at a signal-to-noise ratio 
greater than 5, as in Ref. [6], are masked and inpainted 
using the methods of Ref. [25]. 

The convergence maps are reconstructed from the CMB 
temperature maps as in Ref. [6] using the minimum vari- 
ance quadratic estimator procedure described in Ref. [3]. 
The estimator works as follows: While the unlensed CMB 
is statistically isotropic, any lensing deflection imprints a 
preferred direction into the statistical properties of the 
CMB. This corresponds mathematically to the fact that 
formerly uncorrelated modes of the isotropic unlensed 
CMB temperature field are correlated by lensing, with 
the correlation proportional to the lensing deflection. We 
can hence estimate the lensing convergence by measuring 
the correlation between pairs of Fourier modes using a 
quadratic estimator: 

*(L) = N(L) J d 2 l/(L, l)J(l)r(L - 1). (7) 

Here 1 , L are Fourier space coordinates, N is the normal- 
ization function (which ensures that the estimator is un- 
biased) and / is a weighting such that the signal-to-noise 
ratio on the reconstructed convergence is maximized in the 
case of isotropic noise with full sky coverage (see Ref. [3] 
for the details of these functions). The weighting uses a 
smoothed version of the observed anisotropic noise power. 
To calculate the normalization function, we multiply a 
first-approximation normalization which uses the data 
power spectrum by a small L -dependent correction factor. 
This factor is obtained from Monte Carlo simulations by 
requiring that on average the cross power of the recon- 
structed convergence with the true simulation convergence 
be equal to the true convergence power spectrum. The 
simulated maps used in this Monte Carlo calculation are 
constructed to match ACT data in both signal and noise 
properties, as in Ref. [6]. In our lensing reconstruction, we 
use a wider range of scales in the temperature map than 
in Ref. [6], filtering out modes below € = 500 and above 


€ = 4000. As can be seen in Eq. (7), lensing information at 
a scale L is obtained from two temperature modes sepa- 
rated by L, so that this filtering does not preclude us from 
obtaining small-!? lensing modes. Intuitively, this is 
because we deduce the distribution of large-scale lenses 
from the distortion of small-scale temperature fluctuations. 

We subtract from the reconstructed ACT convergence 
maps a simulated mean field map (k), obtained from 480 
realizations of simulated reconstructed lensing maps. This 
map is nonzero due to correlations induced by window 
functions and noise which, to the convergence estimator, 
appears as a small spurious lensing signal which must be 
subtracted. 

We can thus estimate the CMB lensing quasar cross- 
power spectrum by calculating 

Cf = <Re[(*(l) - </c)(l))*tf(l)]> le€ , (8) 

where the outer average is over all pixels with 
Fourier coordinates 1 which lie within the bandpower 
corresponding to t. 

B. The SDSS quasar maps 

In this work we use the SDSS-XDQSO photometric 
quasar catalog [26,27], extracted from SDSS Data 
Release 8. The analysis used in this catalog separates the 
population of quasars and foreground stars using a proba- 
bilistic model in flux space. This analysis assigns a proba- 
bility of being a quasar to every point source with 
dereddened i-band magnitude between 17.75 and 22.45 
in the SDSS imaging. Though the catalog extends out to 
z > 4, we do not use the highest-redshift sources with 
redshifts z > 3.5. This reduces the shot noise error on our 
measurement, as there are very few such sources. Using 
this catalog, we construct a map of the fractional over- 
density q of quasars across the same 324 square degrees on 
which we perform our lensing reconstruction. We include 
in our quasar maps, with unit weight, all point sources with 
a greater than p = 0.68 probability of being a quasar. As 
the probability distribution is nonuniform, the residual 
level of stellar contamination can be calculated from the 
catalog probabilities to be 6% (we discuss later how this is 
accounted for in our theoretical calculations). A spectro- 
scopic quasar sample at high redshifts (2.2 < z < 3.5) does 
show a contamination fraction which is ~ 15% larger than 
that predicted by the catalog (see Ref. [27]); however, the 
majority of the quasars we consider are at lower redshift, 
where the XDQSO model estimates should be significantly 
more accurate (high-redshift quasar selection is less accu- 
rate, because at z ~ 2. 8 the quasar and stellar loci cross in 
color space). We neglect the error on the stellar contami- 
nation fraction, which is in any case smaller than the 
statistical error (and would only reduce the measured 
signal, not increase it). The area covered by our quasar 
maps contains on average 75 quasars per square degree. 
The redshift distribution of the quasars in our maps is 
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FIG. 1 (color online). The redshift distribution of SDSS qua- 
sars used to construct our maps of fractional quasar overdensity, 
normalized to a unit maximum. The corresponding redshift bins 
are shown with blue filled circles; they are interpolated to give 
the continuous curve used in our theory calculations (blue 
dashed line). For comparison, the red dotted line shows the 
lensing kernel W K (z), again normalized to a unit maximum. 



shown in Fig. 1 , along with the CMB lensing kernel, which 
has a very similar redshift distribution. 

C. The CMB lensing quasar cross-power spectrum 

The cross-power spectrum of the ACT CMB lensing 
maps and the SDSS quasar maps is shown in Fig. 2. 

The error bars on the data points are calculated theoreti- 
cally as proportional to , with an additional factor 

calculated from the number of independent pixels in the 
bin corresponding to each data point. The spectra used in 
this calculation are the full data spectra which include both 
signal and noise (including Poisson noise). Bootstrap error 
estimates from splits of our data are consistent with this 
calculation. For comparison, we also calculate error bars 
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FIG. 2 (color online). The CMB lensing quasar density cross- 
power spectmm, with the data points shown in blue (the covari- 
ance between different data points is negligible). The significance 
of the detection of the cross spectrum is 3.8cr. The green solid line 
is a theory line calculated assuming the fiducial bias amplitude. 
This theory line is reduced by 6% to account for the expected level 
of stellar contamination of the quasar sample. 


using Monte Carlo methods, cross-correlating 480 realiza- 
tions of simulated reconstructed lensing maps (containing 
signal and realistic noise) with the quasar data maps. Both 
methods assume the two maps are uncorrelated; such 
calculations are very good approximations to the true error 
on the cross-correlation because both maps are noisy so 
that C K ( K C qq » (C K f q ) 2 . The error bars obtained from 
Monte Carlo methods are nearly indistinguishable from 
the theory error bars and lead to the same detection sig- 
nificance. (We also verified that replacing the quasar data 
maps with 480 simulated maps with the same number of 
randomly distributed sources leads to similar, though 
slightly smaller, errors.) The Monte Carlo estimates of 
the errors also allow us to calculate the full covariance 
matrix. The off-diagonal elements are negligible compared 
to the diagonal elements; for every bin, the covariance 
between neighboring bins was less than 4% of the bin 
autocorrelation. We thus neglect covariance between dif- 
ferent data points in our analysis. 

Also shown in Fig. 2 is a theoretical calculation of the 
expected cross-power spectrum obtained from Eq. (6). In 
this calculation the matter power spectrum was computed 
using the CAMB software [28]. The nonlinear (HALOFIT, 
Ref. [29]) matter power spectrum was used; however, using 
a linear matter power spectrum instead only slightly 
changed the computed cross spectrum (as most of the 
signal arises from angular scales corresponding to linear 
scales in the matter power spectrum, where the linear and 
nonlinear matter power spectra hardly differ). We use the 
quasar redshift distribution as shown in Fig. 1 in this 
calculation. As the integration kernel is slowly varying, 
the theory curve is insensitive to the binning and interpo- 
lation of this redshift distribution. A fiducial bias model for 
this calculation is obtained by interpolating the central 
measured bias values of Ref. [30] (averaging the values 
obtained with and without the inclusion of negative points 
in the correlation function). The fractional error on these 
central bias values is below 1 0% at low redshifts z < 2, but 
rises to —20% at z — 4-5. The values were obtained from 
measurements of the amplitude of the quasar correlation 
function (which is sensitive to the bias) for a SDSS spec- 
troscopic quasar sample. This fiducial bias model is shown 
in Fig. 3. As the theoretical cross-power spectrum does not 
depend strongly on the detailed form of the bias model, we 
use this measurement as a convenient fiducial template, 
though the spectroscopic catalog used in this measurement 
does not extend to as faint a magnitude as the photometric 
catalog we use to make quasar maps. Despite this, the 
quasar power spectrum predicted by this fiducial bias 
model is consistent with the power spectrum of our quasar 
maps. The calculated theoretical cross-power spectrum is 
reduced by 6% to account for stellar contamination; while 
stars are uncorrelated with lensing, they contribute to the 
average density of sources, and so cause us to calculate a 
fractional quasar overdensity that is 6% too small. 
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FIG. 3 (color online). Blue dashed line: the fiducial quasar bias 
template (interpolated from the data points of Ref. [30]), used in 
the theoretical calculation of the CMB lensing quasar cross-power 
spectrum. Red solid line: the bias amplitude {b/b Rd = 1.02) best 
fit by the measured cross-power spectrum. Red dashed lines: the 
± 1 a error ranges of this amplitude. Note that only one degree of 
freedom is constrained: the overall bias amplitude of an assumed 
redshift evolution. 

The data fit the theory curve (which assumes the fiducial 
bias model) well, with a chi-squared value for this curve of 
Y theory = 13.2 for 10 degrees of freedom. We obtain the 
significance of our detection of the cross-power spectrum 

by calculating jxLn “ theory- where Ynuii is calculated 
for the null line. The significance of the detection is found 
to be 3.8<x. 

IV. A CONSTRAINT ON THE QUASAR BIAS 

We calculate a constraint on the linear bias of quasars 
from the lensing quasar cross-power spectrum. To do so, 
we assume a bias template of the fiducial shape shown in 
Fig. 3, but rescaled by a constant factor for all redshifts. 
We calculate the likelihood as a function of this scaling 
parameter b/b M and plot it in Fig. 4. Our result, b/b tld = 
1.02 ± 0.24, is consistent with the fiducial bias model, i.e., 
a value of unity. (Due to the small negative magnification 
bias, b/b m = 0 does not correspond exactly to the null 
line.) As the redshift distribution is peaked at z = 1.4, 
b/bf jd can be interpreted as approximately parametrizing 
the amplitude of the bias at redshift 1.4. From this inter- 
pretation we obtain a value of the bias at z ~ 1.4 of 
b = 2.5 ± 0.6, which is in good agreement with previous 
measurements from quasar clustering [30]. We can asso- 
ciate this bias with a host halo mass M 200 : using the 
bias model of Ref. [31], we obtain a halo mass of 
log lo (M 2 oo/M 0 ) = 12.9 consistent with previous esti- 
mates [16,30]. We also verify that the cross power calcu- 
lated using only a low (z < 2.2) or high (z >2.2) redshift 
quasar subsample is consistent with the bias given by the 
fiducial model; however, we defer a detailed calculation of 
bias constraints using multiple quasar subsamples to future 
work with a higher signal-to-noise ratio. 
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FIG. 4 (color online). Likelihood as a function of quasar bias 
divided by the fiducial bias, b/b {i& (we assume that the shape of 
the redshift dependence is constant and has the hducial form of 
Fig. 3, and modify the amplitude of the bias function to calculate 
this likelihood.) Interpreting our measurement of b/b ad = 
1.02 ± 0.24 as a bias at z = 1.4 (the peak in the quasar distri- 
bution), we obtain b = 2.5 ± 0.6 at this redshift. 

V. TESTING THE POWER SPECTRUM 
A. Null tests 

We check our result and our pipeline with a number of 
null tests. In a simple first test, we cross-correlate the 
quasar distribution in one part of the sky with the lensing 
convergence in another; as seen in Fig. 5, the results are 
consistent with null as expected, with \ 2 = 6.5 for 10 
degrees of freedom for a fit to null. A more sophisticated 
test is to calculate the cross-correlation of the quasar maps 
with the curl component of the lensing deflection (this 
differs from the convergence reconstructed earlier, which 
is gradientlike). The reconstructed curl map is expected to 
be zero (though it should contain reconstruction noise), and 
hence the cross-correlation with the quasar maps should be 
zero as well. We reconstruct the curl component of the 
estimator as in Ref. [7] [though keeping the normalization 
and filters unchanged from the earlier convergence recon- 
struction, and simply replacing the dot product in /( L, 1) 
of Eq. (7) with a cross product projected onto the 1, X 1, 
direction]. The cross-correlation of the lensing curl com- 
ponent with the quasar maps is also shown in Fig. 5. The 
error bars are calculated from theory as before. As 
expected, this test is consistent with a null result, with a 
(somewhat low) value of x 2 = 3-5 for 10 degrees of 
freedom. 

B. Estimating potential systematic contamination 

We estimate the magnitude of what are expected to be 
the largest contaminants: infrared (IR) sources, SZ clusters 
and galactic cirrus, which contribute flux to the CMB 
temperature maps. (The level of radio source power is 
much smaller, as we can resolve and mask such sources 
down to low flux levels.) Contamination of the cross-power 
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FIG. 5 (color online). Two successful null tests, both consis- 
tent with zero. Upper panel: the cross-power spectrum of quasar 
and lensing maps covering different parts of the sky (permutation 
null test). Lower panel: the cross-power spectrum of the recon- 
structed curl component of the lensing signal with the quasar 
maps (curl null test). 

spectrum is conceivable because the sources of IR and SZ 
signals trace the underlying matter field as the quasars do, 
and galactic cirrus could reduce observed quasar counts by 
extinction. As explained in Ref. [8], any IR or SZ contami- 
nation would appear as a negative bias, as a large IR or SZ 
signal in a certain direction increases local gradients; the 
lensing estimator falsely interprets this as a signature of 
demagnification of the CMB and hence estimates a spu- 
rious underdensity in this direction. 

To obtain an estimate of the contribution of IR and SZ 
contamination to the measured cross-power spectrum, we 
construct simulated quasar maps which are correlated with 
the IR and SZ maps of the simulations of Ref. [32]. The 
quasar maps are constructed by randomly populating all 
halo positions (listed for all halos with masses greater than 
6.8 X 10 12 M g in the catalog supplied with Ref. [32]) with 
quasars. In populating the halos we use a redshift- 
dependent probability of occupation such that the final 
simulated quasar map has the correct redshift distribution. 
Cross-correlating this quasar map with the true conver- 
gence maps of the same simulation, the signal is consistent 


with the theory line of Fig. 2, which confirms that our 
simple simulation has approximately the correct bias (this 
is due to the average mass in the halo catalog being similar 
to the typical halo mass of a linear bias model consistent 
with the fiducial model). 

The level of systematic contamination in our estimator 
can be obtained from these simulated quasar, IR, and SZ 
maps. Keeping the same filtering and normalization as in 
our data, we reconstruct the IR contaminant to the lensing 
signal by replacing the temperature maps in Eq. (7) with the 
simulated IR maps (which we rescaled using an appropriate 
factor as in Ref. [7] to match more recent constraints on IR 
source flux). The cross spectrum of the resulting map with 
the simulated quasar maps gives a negative spurious signal 
which is ~ 7% of the theoretical prediction for the lensing 
quasar cross-correlation. Repeating this analysis with the 
thermal SZ simulations gives a similar negative contami- 
nation of order ~ 5% of theory. (The analysis should over- 
estimate the contamination, as the simulated quasars are 
placed exactly in the centers of the same halos that source 
the SZ and IR signals, neglecting any miscentering effects.) 
Any systematic error in our measurement of the lensing 
quasar cross-power spectrum due to contamination from 
both IR and SZ sources should thus be significantly smaller 
than the size of the statistical error. In addition, the fact that 
these contamination signals are negative means that our 
detection of a positive lensing quasar cross-power spectrum 
cannot be due to such systematics. 

Finally, we estimate the level of contamination from 
galactic cirrus using the dust maps of Ref. [33]. We sub- 
tract a map of the signal at 148 GHz induced by the dust 
(obtained from Ref. [32]) from the ACT temperature data, 
reconstruct lensing, and reestimate the lensing quasar 
cross-power spectrum. We find that the change in the 
cross-power spectrum is very small, of order 3% of the 
theoretical prediction for the lensing quasar cross spec- 
trum. (This is unsurprising, as in our analysis large-scale 
power below € = 500 has been filtered out of the tempera- 
ture maps.) Contamination of the cross power by galactic 
cirrus is thus negligible. 

VI. SUMMARY AND CONCLUSIONS 

In this work we measure the cross-correlation between 
ACT CMB lensing maps and maps of the quasar distribu- 
tion made from the SDSS-XDQSO catalog. We detect the 
cross-power spectrum at 3.8cr significance, directly con- 
firming that quasars trace mass. We check our detection 
with null tests, including a cross-correlation of quasars 
with the reconstructed curl component of lensing, which 
is found to be null as expected. Potential systematic 
contamination is estimated and found to be negligible. 
From our detection we estimate the quasar bias. We mea- 
sure b/b M = 1.02 ± 0.24; interpreting this as a bias at 
z ~ 1.4 (the peak in the quasar distribution), we obtain 
b = 2.5 ± 0.6 at this redshift [which corresponds to a host 
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halo mass of M Q ) = 1 ]. Unlike mea- 

surements from clustering, this lensing measurement 
involves a direct comparison of the quasar distribution 
with the mass distribution, with little modeling required. 

The study of high-redshift mass tracers with CMB lens- 
ing is a new held. In the next few years, the signal-to-noise 
ratio on lensing measurements should improve by an order 
of magnitude with data from experiments such as Planck, 
ACTPol and SPTPol [34-36], ACTPol in particular should 
provide high signal-to-noise-ratio lensing measurements 
which have considerable overlap with SDSS quasar fields. 
Higher signal to noise will allow constraints on quasar 
biases as a function of redshift, luminosity, color or other 
properties and will thus provide a wealth of information on 
the properties of quasars and the halos that host them. More 
precise bias measurements of both quasars and galaxies 
will also allow tests of dark energy properties [37] and 
modified gravity [38]. This work lies at the beginning of an 
exciting research program: the study of astrophysics and 
cosmology with CMB lensing cross-correlations. 
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